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The relativistic hydrodynamical model HYLANDER-C is used to give estimates for single 
inclusive particle momentum spectra in ^/s — 200 GeV/ nucleon Au+Au collisions that will be 
investigated experimentally in the near future. The predictions are based on initial conditions 
that the initial fireball has a longitudinal extension of 1.6 fm and an initial energy density of 
30.8 GeV/ fm? as obtained from a cascade model. For the collision energy considered here, 
different stopping scenarios are explored for the first time. Our calculations give particle 
yields of the order of 10,000 to 20,000 charged particles per event. 
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Forthcoming experiments at the Relativistic Heavy-Ion 
Collider (RHIC) at the Brookhaven National Laboratory 
(BNL) are designed to search for a new state of nuclear 
matter, the quark-gluon plasma (QGP) [1]. The exis- 
tence of a QGP should be reflected in the equation of 
state (EOS) of nuclear matter, which is a relationship 
between extensive variables, e.g., pressure, energy den- 
sity, and baryon density. 

Those theoretical models such as relativistic hydrody- 
namics which are capable of describing relativistic heavy- 
ion collisions and which also allow for an explicit use of an 
EOS have to be based on a thermodynamical approach. 
Several models based on relativistic hydrodynamics exist 
(for recent reviews c/., e.g., Refs. [2,3]); these models as- 
sume the existence of local equilibrium. In general, the 
models solve the hydrodynamical equations - which in 
the case of ideal relativistic fluids are given by the rela- 
tivistic Euler equations [4] - by sophisticated numerical 
techniques. 

There are essentially only three numerical techniques 
that are employed in the study of heavy- ion reactions [2] : 
particle methods, continuum methods, and predictor- 
corrector methods. In our attempt to make predictions 
for single inclusive particle momentum spectra of various 
hadron species emerging from ^/s = 200 GeV/nucleon 
Au-I-Au colhsions, we shall use a predictor-corrector 
model, which has been successful in providing a consis- 
tent and coherent description of the single and double 
inclusive cross sections of mesons and baryons for several 
types of flxed target heavy-ion collisions at CERN/SPS 
beam energies near 200 A-GeV: HYLANDER-C [5], 
which is the upgraded successor of HYLANDER [6] . 



*E. Mail: schIei@LANL.gov 
tE. Mail: dds@LANL.gov 



Predictor-corrector type fluid models begin the calcula- 
tions with the matter already highly compressed and the 
numerical code then follows the expansion of the fluid as 
governed by a particular EOS. At relativistic energies, it 
is important to employ a covariant description of freeze- 
out, which allows the calculation of the various particle 
spectra [7] . In earlier work, the initial state of the hot 
and dense zone of nuclear matter - the fireball - was 
parametrized, and was adjusted until a match with ob- 
servables was established [8]. HYLANDER-C, e.g., can 
use initial conditions ranging between the extremes de- 
fined by the Landau [9] (i.e., complete stopping) and the 
Bjorken [10] (i.e., minimum stopping through a scaling 
ansatz) initial conditions. 

Conversely, one could also take information about the 
initial state of the fireball from other sources, e.g., a cas- 
cade code, as input and then allow the system to evolve 
according to an assumed EOS, and finally particle mo- 
mentum spectra will be calculated. As there are no data 
yet available for — 200 GeV/nucleon Au-|-Au col- 
lisions, we shall discuss in this paper results for a rel- 
ativistic heavy-ion collision at these high energies using 
HYLANDER-C, and the sensitivity of observables to dif- 
ferent assumptions regarding different initial conditions. 
For the following we have to specify an equation of state 
and the initial conditions. In addition we choose for 
our freeze-out condition a fixed freeze-out energy den- 
sity, e/, and we assume that the freeze-out occurs for all 
particle species at the same fixed value of e/. The cal- 
culations using HYLANDER and HYLANDER-C which 
reproduced data of the CERN NA35, NA44 and NA49 
collaborations [11] - [13] were obtained while using an 
EOS with a phase transition to a QGP at a critical tem- 
perature Tc = 200 MeV {cf. Refs. [14,15], and Refs. 
therein). The EOS used in this work does not depend 
on the baryon density (which is a valid assumption for 
the collision energy regime considered here) , and thus the 
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freeze-out energy density translates into a fixed freeze-out 
temperature Tf. The choice for the freeze-out tempera- 
ture has been Tf = 139MeV in the former calculations; 
we shall use the same value here. This is consistent with 
our attempt to be conservative in determining possible 
results for the Au+Au collisions by choosing exactly the 
same conditions for the fireball expansion (through the 
same choice of the EOS) and the same conditions for the 
freeze-out (through the same choice for e/ and Tf) that 
were used in the successful description of the CERN/SPS 
data. 

Because the appropriate initial conditions for the forma- 
tion of the Au+Au fireball are unknown, we shall adopt 
values for the hydrodynamical input from K. Geiger's 
parton cascade calculations [16-18] (other cascade mod- 
els [19,20] have been employed for predictions of the 
heavy-ion collision discussed here with results that are 
similar to K. Geiger's results). In Ref. [17] it was re- 
ported that after an equilibration proper time, Teq = 1.8 
fm/c, the initial fireball originating from a ^/s = 200 
GeV/nucleon Au-I-Au collision was formed with an ini- 
tial longitudinal extension (of the Landau volume [21]), 
A = 1.6 fm, and an upper value for the initial energy 
density, eA = 30.8 GeV/fm^. 

We use here for the initial distributions the initial condi- 
tion scenario which has been described in Refs. [21,22]. 
Specifically, our model uses the five initial parameters, 
Kl, a, j/a, Vm and cr, which are the relative fraction 
of thermal energy in the central fireball, the longitudinal 
extension of the fireball, the fluid's rapidity at the edge 
of the central fireball, the fluid's rapidity at the maxi- 
mum of the initial rapidity distribution of the baryons, 
and the width of the initial baryonic rapidity distribu- 
tion, respectively. Furthermore, it is assumed that an 
initial transverse fluid velocity component is absent, and 
the initial longitudinal distributions (with respect to the 
beam axis) for energy density, e, and baryon density, ub- 
arc smeared out with a Woods-Saxon paramctrization in 
the transverse direction, r± {cf. Refs. [22,23]). 
In the following, we keep A and Kl fixed to the values 
provided by K. Geiger's results. We stress that an initial 
energy density, eA = 30.8 GeV/ fw? in the Landau vol- 
ume represents only 11.5% of the total available energy. 
The parameter a we keep fixed also without loss of gener- 
ality. In particular, it is our intent to discuss six possible 
stopping scenarios {cf. Table I, scenarios I - VI). We can 
accomplish this by varying the only two parameters left in 
our model for the initial conditions, j/a and Dm- A larger 
value for the absolute value of rapidity, t/A, at 2: = ±A/2 
results in a larger initial rapidity field, y{z), which can be 
interpreted in having less stopping. Conversely, initially 
concentrating as much baryonic matter as possible in the 
central region of the initial fireball as energy conserva- 
tion permits us, i.e., by minimizing y™, results in having 
a much larger stopping compared to the cases where the 
baryonic matter lies further outside. For the following 



discussion, wc have chosen three scenarios where as much 
baryonic matter as possible is concentrated in the central 
region of the initial fireball (scenarios I, III, V), and three 
scenarios where the baryonic matter lies further outside 
(scenarios II, IV, VI). 

In Fig. 1 we show the initial distributions of energy den- 
sity, e, baryon density, n^, and the fluid rapidity, yp, 
normalized to their maximum values, e™'^^, B'^^^, and 
Hem, and plotted against the longitudinal coordinate z, 
for each of the six initial conditions. The absolute values 
of the initial fluid rapidity, j/Aj &t z = ±A/2 decrease 
from the top to the bottom of Fig. 1. Therefore, we have 
a greater amount of stopping as one goes from the top to 
the bottom of Fig. 1. Comparing the left column of Fig. 
1 (scenarios I, III, V) with the right column (scenarios 
II, IV, VI), we recognize that we have greater stopping 
when going from the right to the left in the figure. Table I 
gives also the relative fractions, , of baryons initially 
concentrated in the initial fireball volume of length A. 
Thus, scenario V gives the highest initial stopping, and 
scenario II gives the smallest initial stopping. 
Inspecting Fig. 1 further, we also see that the initial 
conditions I and II are closer to a Bjorken-type initial 
condition, and scenarios V and VI are closer to a Landau- 
type initial condition. 

After fixing the various initial conditions, we allow for 
expansion and cooling of the relativistic fluids by solving 
the relativistic Euler equations. As mentioned above, 
freeze-out occurs for all particle species at the same fixed 
value of Tf = 139 MeV. In Fig. 2 we display for 
each scenario the isotherms of the expanding systems 
at the radial coordinate, r± = 0. In particular, the 
outer isotherms represent the freeze-out hypersurfaces at 
rj_ = 0. The total lifetimes, tmax, of the various systems 
range from 37.8 fm/c to 49.7 fm/c, and the total life- 
times of the QGP, tgop, range from 16.9 fm/c to 23.7 
fm/c. We note that tqciP is not correlated with tmax be- 
cause of the complex behaviour of the numerically treated 
multi-dimensional relativistic fluids. Furthermore, look- 
ing at Fig. 2, we can see that the freeze-out hypersur- 
faces are to a large extent hyperbola-like shaped, which 
we especially would expect from a purely Bjorken type 
expansion scenario [10]. 

Following the formalism which is outlined in refs. [21,22], 
we have calculated single inclusive particle momentum 
spectra of pions, kaons, and protons for ^/s = 200 
GeV/nucleon Au-|-Au collisions. We have taken into ac- 
count resonance decay contributions up to the third gen- 
eration of particle production. Figs. 3 and 4 show our 
results for rapidity and transverse momentum spectra, re- 
spectively, for the six initial conditions. It is remarkable 
that a wide range of shapes is possible for the rapidity 
spectra. We observe an almost Gaussian distribution for 
scenario V, a flat rapidity plateau for scenario III, and 
saddle-shaped distributions for all the other four scenar- 
ios. In agreement with expectation, the rapidity spec- 
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tra are narrower and the transverse momentum spectra 
show a higher mean transverse momentum, {k±), (i.e., a 
smaller slope,) for initially higher stopping. Contrary to 
the common belief [10], Bjorken type expansion scenarios 
do not necessarily lead to flat rapidity plateaus. This is 
because of finite size effects of the fireballs. 
In addition, we list in Table II the absolute particle yields 
for some selected particle species. The numbers were 
obtained by integrating the corresponding rapidity and 
transverse momentum spectra, respectively. We stress 
that the particle yields in Table II are obtained while 
assuming that ca = 30.8 GeV/fm^, i.e., = 11.5%. 
If one would choose pure Landau initial conditions (i.e., 
complete stopping) one would obtain approximately 10 
times that many particles. 

To summarize, we have used the relativistic hydrodynam- 
ical model HYLANDER-C to give estimates for single 
inclusive particle momentum spectra of various hadrons 
for ^/s = 200 GeV/nucleon Au+Au collisions. Our pre- 
dictions are based on the assumptions that the initial 
fireball has a longitudinal extension of 1.6 fm and an 
initial energy density of 30.8 GeV/ fm^ , values consistent 
with results obtained using the cascade code of Geiger. 
Within different stopping scenarios, we have investigated 
the freeze-out hypersurfaces and the corresponding fi- 
nal particle distributions. Our calculations give parti- 
cle yields of the order of 10,000 to 20,000 charged par- 
ticles per event. We have also found that the shapes of 
the rapidity spectra are sensitive to the assumed stop- 
ping power. Our quantitative estimates are based on 
assumptions for the initial conditions obtained from cas- 
cade models, the EOS, and the proper freeze-out condi- 
tions for the heretofore experimentally unexplored regime 
of ^/s = 200 GeV/nucleon collisions. Information con- 
cerning each of these must be gained with additional the- 
oretical work as well as data when they become available. 
We iU'Q grateful for many instructive discussions with Dr. 
Klaus Kinder-Geiger. This work has been supported by 
the U.S. Department of Energy. 
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TABLE I. Properties of the initial fireballs. 
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11 


III 


IV 


V 


VI 




Initial parameters 








Kl 


0.115 


0.115 


0.115 


0.115 


0.115 


0.115 


A [fm] 


1.6 


1.6 


1.6 


1.6 


1.6 


1.6 


J/A 


0.9 


0.9 


0.6 


0.6 


0.3 


0.3 


Vm 


1.52 


2.00 


1.24 


1.75 


0.88 


1.50 


a 


0.4 


0.4 


0.4 


0.4 


0.4 


0.4 


Tf [MeV] 


139 


139 


139 


139 


139 


139 






Output 








ycm 


5.36 


5.36 


5.36 


5.36 


5.36 


5.36 


eA [GeV/fm^] 


30.8 


30.8 


30.8 


30.8 


30.8 


30.8 


^max [GeV/fm^] 


40.7 


30.8 


45.8 


30.8 


41.8 


30.8 




2.23 


2.09 


1.52 


1.44 


0.78 


0.74 


JriB 


0.062 


0.003 


0.055 


0.002 


0.072 


0.001 


tmax [fm/c] 


49.7 


51.7 


45.6 


49.0 


37.8 


43.5 


tQGP [frn/c] 


23.7 


16.9 


23.5 


20.3 


21.2 


20.6 


TABLE II. Particle yields of pions, kaons, and proto 


ns. 


I 


II 


III 


IV 


V 


VI 


5320 


4654 


5802 


5501 


6816 


5741 


N^- 5420 


4748 


5903 


5609 


6932 


5840 


Nji+ 770 


674 


834 


798 


972 


822 


Nk- 709 


616 


776 


735 


915 


766 


Np 328 


305 


336 


352 


364 


333 
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FIG. 1. Initial distributions of energy density, e (solid 
lines), and baryon density, ub (dotted lines), as well as the 
fluid rapidity, yp (dashed lines), normalized to their maxi- 
mum values (cf. Table I) and plotted against the longitudinal 
coordinate z. 




Rapidity 

FIG. 3. Rapidity spectra of negative pions, 7r~ (solid lines) , 
negative kaons, K~ (dotted lines), and protons, p (dashed 
lines). In each plot the rapidity spectra of K~ and p are 
enhanced by a factor of 5 for better visibility. 
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FIG. 2. Thermal evolution of the Au+Au fireballs. In each 
case the outer contours are isotherms for T = 140 MeV , and 
each successively smaller contour represents an increase of 20 
MeV in the temperature. 
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FIG. 4. Transverse momentum spectra of negative pions, 
7r~ (solid lines), negative kaons, K~ (dotted lines), and pro- 
tons, p (dashed lines). 
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